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O ■ ABSTRACT 

^ ' We present very deep and accurate photometry of the open cluster M35. We have 

observed this association in the Cousins R,I filters, together with the Johnson V filter. 
^ ' We have covered a region of 27.5x27.5 square arcmin, equivalent to a fifth of the total 

■ area of the cluster. The data range from Ic=12.5 to 23.5 magnitudes, and the color 

intervals are 0.4<(V-I)c<3.0, 0.5<(R-I)c<2.5. Roughly, these values span from 1.6 
down to the substellar limit, in the case of cluster members. By using the location 
of the stars on color- magnitude and color-color diagrams, we have selected candidate 
members of this cluster. We have merged our sample with previously published data and 
obtained a color-magnitude diagram for the complete stellar population of the cluster. 
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covering the spectral range early B - mid M. Based on the distribution of field and 

cluster stars in color-magnitude and color-color diagrams, we estimate that two thirds 
of these candidates are likely to be true members of M35. These stars approximately 
double the number of stars identified as candidate members of this cluster (~2700). We 
provide the photometry and accurate positions of these stars. The deep photometry has 
allowed us to study the mass segregation within the cluster, the Luminosity Function 
and Mass Function. We show that in the magnitude range 13 < Ic < 22 there is a 
reduced mass segregation, in opposition to what happens to higher mass stars, where 
the mass segregation is stronger. The Luminosity Function behaves essentially as the 
one characteristic of the Pleiades, presenting a peak at Ic ~19 magnitudes (Mj ~9). 
Combining our photometry with previous data corresponding to more massive stars, we 
find that the Mass Function increases monotonically, when plotted in a log-log form, 
until it reaches ~0.8 M0 (q:=2.59). It remains shallower for less massive stars (a=0.81 
for 0.8-0.2 Mq), whereas a decrease ins observed for stars close to the substellar regime. 
These different behaviors suggest that at least three mechanisms play a role in the 
formation of stellar and substellar objects. The total mass of the cluster is ^^1600 Mq 
in the area covered by this study. 

Subject headings: open clusters and associations: M 35, NGC 2168; photometry; late 
spectral type stars; stellar mass function 

1. Introduction 

Open clusters serve as Rosetta stones to study stellar properties because they provide homo- 
geneous samples of known age, metallicity and distance. Therefore, open clusters can be used to 
determine empirically the evolution of different stellar properties such as the activity, rotation, 
chemical abundances, mass loss, et cetera. For low mass stars, however, only a handful of open 
clusters have been studied systematically, essentially due to the fact that most clusters are far away, 
with distance moduli ranging from 6 to 12 magnitudes. With the recent improvements in detectors 
and with the availability of large ground-based telescopes, it is possible to obtain high quality data 
for a larger number of clusters. 

M35 (NGC 2168) is a very interesting open cluster for several reasons. It is moderately nearby, 
with (m-M)o=9.7 (Vidal 1973). It is very well populated, containing several thousand members 
with a total dynamical mass estimated between 1600 and 3200 Mq (Leonard & Merritt 1989). It is 
spread over an area of ~lxl square degrees, making it a suitable target for multifiber spectroscopy. 
The contamination by field stars is not overwhelming, as happens in some other clusters of similar 
characteristics. The reddening is also moderate -E(B-V)=0. 17-0.225. Perhaps one of the most 
interesting properties is its youth, since it is said to be a cluster coeval with the Pleiades, based 
on the fact that their upper-main sequence turn-off locations are said to be at the same position in 
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color-magnitude diagrams (Vidal 1973). 

Until very recently, only the upper-main sequence (MS) of M35 had been studied with some 
detail. Photometric (UBV filters) and/or proper motion searches were carried out by Cuffey (1938), 
Hoag et al. (1961), Cudworth (1971), Vidal (1973), McNamara & Sekiguchi (1986a). They reached 
barely V=15, about 1.2 Mq, in the case of M35 members. At the present time, different groups have 
conducted several multi-band CCD photometry which reach much deeper members, such as Sung 
& Bessell (1999=SB1999), von Hippel et al. (2000) and Sarrazine et al. (2000). The last works 
present data in the UBVRI filters, reaching beyond V=19, whereas Sung k, Bessell (1999) used the 
UBVI passbands only, reaching one magnitude fainter. Sung & Bessell (1999) were able to estimate 
the interstellar reddening, the distance and the age of the cluster, deriving E(B-V)=0.255ib0.024, 
(m-M)o=9.60±0.10 and r=200lioo Myr. Previously, Sung & Lee (1992) derived a distance modulus 
(m-M)o=9.3, an age r=85 Myr, and a differential internal extinction E(B-V)=0.26-0.44. A different 
estimate of the distance and reddening have been published by von Hippel et al. (2000), E(B- 
V)=0.3, (m-M)y=10.25. These authors also estimated an age by fitting isochrones, yielding r=100 
Myr. However, they inferred that M35 is older than the Pleiades by comparing in a differential way 
members of both clusters. Based on an extension of the data presented by von Hippel et al. (2000), 
and assuming a metallicity of [Fe/H]=-0.21 (Barrado y Navascues et al. 2000a, b), Sarrazine et al. 
(2000) derived (m-M)=10.16±0.10, E(B-V)=0.198±0.008 and r=160±40 Myr. This age and the 
Sung &; Bessell (1999) value are significantly older than the canonical upper-main sequence age 
of the Pleiades, 70-100 Myr (Mermilliod 1981; Meynet, Mermilliod, & Maeder 1993). These older 
ages are in agreement with our own findings concerning the distribution of rotational velocities 
and lithium abundances (Barrado y Navascues et al. 2000a, 2000b). Based on all the available 
information, they assumed an age of r=175 Myr. Note that Barrado y Navascues ct al. (2000b) 
have derived the first metallicity estimate for M35 based on high resolution spectroscopy, with 
the resulting value of [Fc/H] = -0.21ib0.10, lower than the value corresponding to the Pleiades, 
[Fe/H]=-|-0.01, estimated in the same fashion by these studies. The Pleiades cluster is one of 
the cornerstones of astrophysics, and the comparison of members of both clusters can provide 
additional information for understanding their properties. Although several stellar characteristics 
are thought to depend essentially on age (rotation, coronal and chromospheric activity, et cetera), 
this assumption has been recently challenged by Randich k. Schmitt (1995), comparing the X-ray 
luminosities of Praesepe and the Hyades, two coeval clusters (see also the discussion in Barrado y 
Navascues et al. 1997, 1998a). Extensive and detailed studies of several clusters of similar age can 
help to understand this crucial problem, as well to clarify the evolution of the stellar properties 
when comparing clusters of different ages. 

We have embarked on an ambitious study of the open cluster M35, combining high resolu- 
tion spectroscopy to determine lithium abundances, radial and rotational velocities (Barrado y 
Navascues et al. 2000a, 2000b) with deep imaging to extend the membership list to nearly the 
hydrogen burning mass limit (or beyond) and estimate the cluster Mass Function (Barrado y 
Navascues et al. 1998b). The current paper is devoted to the presentation of our deep imaging 
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data, to study the complete MF, the mass segregation within the cluster and to estimate the total 
mass of M35. 



2. Observations 

We have carried out photometric observations of the young cluster M35 in two different cam- 
paigns. The first run took place in December 11th, 1995, at the KPNO^ 4m telescope. In this 
case, we used the prime focus, a CCD camera, with a Tek 2048x2048 pixels, yielding a field of 
23.2x23.2 square arcmin. For the second run we used the Canada-France-Hawaii Telescope (De- 
cember 13th, 1996) and the mosaic camera (8000x8000 pixels), which provided a field of 27.5x27.5 
square arcmin. 

The observations were processed using standard procedures (bias subtraction, flat-field cor- 
rection, etc). We used the IRAF^ package in the case of the KPNO data. Since our aim was to 
obtain photometry of cool stars over a very large dynamical range, we collected two sets of images 
of different exposure times during the KPNO run, namely 60 and 600 sec for the V filter, and 10 
and 120 sec for the Cousins Ic band (Cousins 1978). Therefore, we reached (V,Ic)iimit — 

(20.2,18.2) 

for the shallow set and (V,Ic)iimit=(21-l,19.5) for the deep exposures. The completeness limits are 
Vcompiete=19-0 ^nd Vcompiete=19-75, respectively. In the case of the CFHT survey, we collected 
data in the magnitude interval 17.5<Ic<23.5, 0.5 < (R-I)c < 2.5. The completeness limit was 
Ic, complete =22.1 magnitudes in this case. Errors can be estimated as 0.05 magnitudes for each filter 
at the bright end, degrading to 0.15 magnitudes for the faintest objects in each dataset. 

For the KPNO observations, the calibration of the data was performed by observing standard 
stars in Landolt (1992) -SA 92, 94, 95, 98, 100, 112, and 113- several times during the observing 
run. We derived our own extinction coefficients, yielding Cv=0.192±0.006 and Ci=0.110±0.003. 
The calibrated data was obtained by using the expressions: 

V = v- 3.583 - 0.192 x airmass - 0.046 y. {v - i) (1) 

I = i- 4.276 - 0.110 x airmass + 0.012 x {v - i) (2) 

where V,I denote the calibrated magnitudes, whereas v,i refer to the instrumental values. 

Wc have compared our final calibrated data with the two previous studies which arc deep 
enough. The comparison between our data and the values from Sung & Bessell (1999) indicates 
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that the I data are essentially identical, whereas there is a difference in the data corresponding to 
the V, in the sense V(SB1999)=V(this work)+0.03 magnitudes. Of this difference, 0.01 magnitudes 
can be attributed to the difference between the Landolt (1992) photometry (the values used for 
the calibration), and the SAAO photometry, which defines the Cousins system. Regardless of this 
systematic difference, our final results (membership selection and Mass Function analysis) are not 
affected by it. 

The CFHT observations were carried out and analyzed in the same fashion as our observations 
of the Pleiades and Praesepe (Bouvier et al. 1998). In all cases the data were recorded under 
photometric conditions and the calibration was performed independently from each other using 
standard stars listed in Landolt (1992). Several very red standard stars were observed to ensure the 
accurate transformation from the instrumental system into the Cousins system. Visual inspection 
of the stars in common to the two set of data (KPNO and CFHT), as listed in Table 2, indicates 
that the calibrated magnitudes are in excellent agreement. 

In the case of the KPNO survey, coordinates were derived for all the detections using the 
GSC catalog. Positional accuracies are better than 0.5 arcsec. An empirical confirmation for the 
accuracy of the coordinates was provided by the good quality of the spectra we obtained using 
our coordinates for a WIYN/HYDRA multi-object spectroscopy run (Barrado y Navascues et al. 
2000a, 2000b). 

3. The Selection of Members 

3.1. A Main Sequence for Cool Stars 

In order to select possible M35 cluster members based on their location in a color-magnitude 
diagram (CMD), we need to have a reliable main sequence locus for low mass stars. Theoretical 
isochrones (e.g. Baraffe et al. 1998) do not yet do this optimally for optical colors due to remaining 

problems with cool dwarf model atmospheres and inadequate detailed line opacities used to compute 
fluxes. Therefore, we must rely on an empirical main sequence. We have constructed one using 
field stars from Leggett (1992) which fulfill two conditions: they have accurate parallaxes (11), in 
the sense (T(n)/n<0.10, and they were not listed as "AB" in that paper (that is, they are not 
spectroscopic binaries or close visual binaries). The selected data come from Luyten (1972), Giclas 
et al. (1971), Gliese (1969) and Glicsc & Jahrciss (1979) -see reference therein for the original 
sources of the photometric data-. Figure 2 of Bouvier et al. (1998) and Figure 1 of Stauffer et 
al. (1999) show this empirical MS when applied to the Pleiades and Alpha Per open clusters, 
respectively, in the [Ic,(R-I)c] plane. Figure 1 of Barrado y Navascues et al. (1999) displays the 
comparison between the empirical MS and field stars from the Leggett (1992) in the [V,(V-I)c] 
plane. As can be seen, the lower envelope is an appropriate MS. 

For the current paper, we need to extend this MS to bluer (higher mass) stars. We have done 
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this first by tracing a lower envelope to M35 members identified by Sung and Bessell (1999) in the 
Ic magnitude range 8.5 to 14.0. Second, we used accurate radial velocities for 78 of our photometric 
candidates down to Ic=16, to determine whether they are indeed members (Barrado y Navascues 
et al. 2000b, shown as solid circles). The loci of those bona fide members were used to extend the 
Sung & Bessell (1999) data and to connect it with the empirical isochrone derived from Leggett's 
data. The final result can be seen in Figure la and Figure lb. Note that for the lowest masses, we 
would expect M35 members to be above (to the right) of this curve because the M35 members will 
still be contracting to the main sequence, assuming the canonical age of the Pleiades. However, the 
low mctallicity of M35 -[Fe/H]=-0.21-, produces a shift of the location of its members compared 
similar stars of solar composition and as discussed in Section 1, there are evidences that M35 is 
older than the Pleiades. At higher masses, our MS curve should serve as an approximate lower 
envelope to the distribution of M35 members (e.g. photometric binaries should still scatter up to 
0.75 mag above the curve). Our empirical MS can be found in Table 1. 



3.2. Color- magnitude and color-color diagrams 
3.2.1. KPNO Data 

Figures la and lb depict the photometric Vic data collected at KPNO. The first panel shows 
the CMD for the shallow exposure, whereas panel b contains the values for the deep one. In both 
cases, limiting and completeness magnitudes are indicated as solid and dashed lines. Our empirical 
MS sequence is also shown as a solid line. We used a distance modulus of (m-M)o=9.60, a color 
excess of E(V-I)c=0.321 and an interstellar absorption of Av=0.765 (Ry=3.0), which corresponds 
to the reddening estimated by Sung &; Bessell (1999) of E(B-V)=0.255. The probable members of 
M35 (Barrado y Navascues et al. 2000b), observed spectroscopically (<RV>=-8 km/s), are shown 
as solid circles. 

We have selected as initial possible candidate members of M35 all those stars which have 
photometry in a strip described by the empirical MS and the MS minus 0.75 magnitudes (to allow 
for binarity) and ±0.1 magnitudes, to allow for photometric errors and errors in the distance and 
the reddening. Hence, the strip is 0.95 magnitudes wide. Another example of photometric selection 
of candidate members of a cluster can be found in Barrado &: Byrne (1995), in the case NGC5460, 
and Barrado y Navascues et al. (2000c) for IC2391. The selection of M35 candidate members was 
stopped at V=21.0 magnitudes. Therefore, we have covered the spectral range F0-K7 for cluster 
members. Note that despite the lack of separation between field stars and the locus of the M35 
main sequence, the diagrams show a relative concentration of objects corresponding to M35, very 
clear in Figure lb in the range 1.5<(V-I)c<2.0. 
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3.2.2. CFHT Data 

Our initial optical photometric survey for members of M35 was extended by taking longer 
exposures at the Canadian-Prench-Hawaii telescope. In this case, we obtained 2 exposures in the 
Cousins R and I filters. Figure Ic displays the Ic versus (R-I)c color-magnitude diagram. The 
objects detected by this survey are plotted as dots. Wc have added the empirical MS for field stars 
(see section 3.1), shifted according to the reddening -E(R-I)c=0.178- and the distance modulus 
-(m-M)i = 10.044- of M35 (see previous subsection). Note that the MS really traces the relative 
concentration of very cool stars in M35, although there is a shift of ~0.1 magnitudes, which can 
be attributed to the low metallicity of the cluster. Note that the lower, redder end of the CMD 
corresponds to objects of M5.5 spectral type. The location of the M35 low mass candidate members, 
below our empirical MS, seems to indicate that they are not pre-MS objects, and therefore, they 
are older than ~100 Myr, in agreement with the age by Sung &: Bessell (1999) and Sarrazine et al. 
(2000) and our previous findings (Barrado y Navascues et al. 2000b). 

Figure Ic clearly shows two groups of stars: A significant fraction of the stars are on the locus 
of M35. Most of the detected stars lie between 0.5 < (R-I)c < 1-7, following a distribution parallel 
to the M35 MS, but 4 magnitudes fainter. However, there is not a clear gap, empty of stars, between 
the former (M35 candidate members) and the latter groups (field stars). A smooth transition is 
present, indicating that any list of candidates of M35 presents contamination by field stars. In any 
case, the diagram allows us to perform a tentative identification of candidates for membership to 
the cluster. We have identified candidate members of M35 here in the same way that we did for 
the KPNO data, except in this case allowing also for the fact that the faintest M35 members will 
still be contracting to the MS, and for uncertainties in our empirical MS, as well as in the cluster 
age, distance and reddening. We have cross-correlated the list of stars selected from the Vic plane 
with the list coming from the RIc survey, merging both sets of data. Note that the intersecting 
magnitude range is only one magnitude wide, from 17.5 to 18.5 in Ic. 

3.2.3. The (R-I)^-(V-I)^ Plane 

Our two surveys (CFHT and KPNO) overlap over a range of one magnitude. Therefore, we 
have merged both datasets, generating a subsample of data with measurements in three passbands. 
Figure 2 shows (R-I)c against (V-I)c for this subgroup. We have included in the figure the average 
behavior for field stars, after Leggett (1992) and Bessell (1991), which are shown as solid and 
dashed lines, respectively. Very late spectral type members of the young clusters IC2391 (Barrado 
y Navascues et al. 2000c) are also included as a comparison group. Based on this diagram, we have 
rejected several stars as members of the cluster, due to their location (cross symbols). The solid 
dots represents the photometric members of the clusters. 
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3.2.4- A List of Members 

Our list of members contains 1945 stars, with 275 stars previously observed (essentially by 
McNamara & Sekiguchi 1986a and Sung &; Bessell 1999). Therefore, we present a list with 1670 
new possible photometric members of the cluster, doubling the number of candidate members of the 
cluster. Including all previously known members of M35, there are in total about 2710 identified 
candidate members. As a comparison, the Pleiades has 1194 known members (Pinfield et al. 1998). 

Table 2 lists our identifying number for each object and the photometry in columns #1 to #6. 
Positions are listed in columns ^7 and ^8. Finally, cross identification with WEBDA database 
(Mermilliod 1996), McNamara & Sekiguchi (1986a) and Sung &: Bcsscll (1999), and membership 
probability, derived from the proper motion, are included from column #9 to column #12. Note 
that Mermilliod's identification number in his WEBDA database coincide with the identification 
number of Cuffey (1938) for the first 778 stars. Stars in the M35 field of view in WEBDA database 
are identified with numbers from 1 to 4219 (the database includes 1566 entries, some of them are in 
fact no members) . We have decided to identify our M35 candidate members with numbers starting 
at 5001. 

Our final membership assignment is also listed in Table 2. Those stars with radial velocity 
and/or with proper motions in agreement with the cluster average are flagged with "Y+" in column 
#13. The same criteria were used to reject several candidates and they have the "N+" flag. 
Candidate members with positions in the color-color diagram which do not correspond to the 
cluster locus appear as "N", whereas those stars which only have V,Ic and/or Rc,Ic data are listed 
as "Y". 

Final color-magnitude diagrams, in absolute Ic magnitude and dereddened colors, are presented 
in Figure 3a, Figure 3b and Figure 3c. In all three cases, our empirical main sequence is displayed 
as a solid line. In Figure 3a, data from the literature (Sung & Bessell 1999) appear as asterisk 
symbols (membership probability, from the proper motions, larger than 0.80, "Y+"). Our bona 
fide members, selected from their radial velocity (Barrado y Navascues et al. 2000b), are shown as 
solid circles ("Y+"). Other photometric members are displayed as open circles ("Y"), whereas non- 
members appear as crosses ( "N" , "N+" ) . Figure 3b displays the photometric candidate members 
from the CFHT survey -Ic versus (R-I)c- A complete color-magnitude diagram for all bona fide and 
candidate members of M35 is shown in Figure 3c, stretching from M(Ic)=-1.5 up to M(Ic)=13.0, 
and covering from early B down to M5.5 spectral types. It includes our own candidate members 
and those candidate members from Sung & Bessell (1999) which are probable members, based on 
their proper motions, as derived from McNamara & Sekiguchi (1986a). Note that in order to show 
our CFHT data at the red, faint end, we derived (V-I)c colors from (R-I)c and our empirical MS. 
The star HD41996, which is a likely member based on its proper motion and has already evolved 
off the MS, is out of the limits of the figure. 
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3.2.5. Brown Dwarfs in the cluster 

Our fainter candidates reach Ic~23.0, which corresponds to M(Ic)~13.0 magnitudes. Following 
Stauffer, Schultz &; Kirkpatrick (1998) and the theoretical models of Baraffe (1998, priv. commu- 
nication), the substellar limit should be, for an age of 125-200 Myr, at Ic(BD)~22. 1-22.6, and 
the lithium depletion boundary (sec Stauffer & Barrado y Navascucs 2000) at Ic(LDB)~22.3-23.2 
magnitudes. Therefore, we have reached the substellar domain in M35, and found a sample of the 
farthest brown dwarf (BD) candidates discovered so far. In total, we have 65 objects with masses 
in the range 0.075-0.055 M©, based on a 100 Myr isochrone from D'Antona & Mazzitelli (1994). 
Note that if the masses of cluster brown dwarfs sum up about a few percent of the total mass of the 
association, it would mean that M35 has a minimum of 600-1000 substellar objects. For instance, 
brown dwarf searches in the Pleiades have estimate that the total contribution of BD to the mass 
of the cluster can be 5% (Bouvier et al. 1998) or 3% (Hodgkin & Jameson 2000). 

3.3. Contamination by field stars 

The relatively smooth transition between the concentration of field stars, fainter and/or bluer 
than M35, and this cluster, indicates that our list of member candidates is contaminated by field 
stars. However, note that in the case of the V,(V-I)c diagram, the gap between the M35 main 
sequence and the bulk of the field stars increase with the magnitude: the separation in color is 
smaller at the upper part of the diagram and, therefore, the number of spurious members should 
be larger for stars of spectral type G and early K. Moreover, the separation in color between the 
field stars and the cluster locus is larger in the case of the I,(R-I)c diagram, a fact which indicates 
that it is easier to estimate the contamination by field stars in this case than from the V,(V-I)c 
diagram. 

In the case of the brighter stars presented in this work, which come from the V,(V-I)c CMD, we 
have estimated the contamination by field stars using two different methods: (i) Some of the stars 
have been discovered previously, and proper motions and their associated membership probability 
are listed by the authors (Cudworth 1971; McNamara Sz Sekiguchi 1986a). They have magnitudes 
in the range 13.7<V<14.9, which correspond to F2-G2 spectral types. We have compared the 
number of stars present in our survey having probability less than 50% with the total number, and 
assumed that this ratio is the contamination rate. The result is 37%. If a more demanding criterion 
is selected (such as membership probability less than 80%), the pollution rate does not change much 
(53%). (ii) We have previously derived radial velocities for a subsample of these stars (Barrado 
y Navascues et al. 2000b). They have magnitudes in the range 14.5<V<17.5, corresponding to 
F9-K5 spectral types. We have used the radial velocity as a membership criterion. The pollution 
rate derived from the comparison between the bona fide members and non-members is 30%. These 
numbers confirm our initial conclusion that contamination is stronger in the upper Main Sequence 
(for F and G stars) than for K members. 
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in the case of the deep survey from CFHT, we have divided the Ic,(R-I)c color-magnitude 

diagram in boxes of size 0.5 magnitude by 0.1 magnitude in Ic and (R-I)c, respectively. We have 
constructed diagrams representing the number of stars in each i box -Nj^*- against (R-I)c for a 
given Ic interval. Figure 4 shows one such diagram, corresponding to the interval 18.0 < Ic < 18.5 
magnitudes. The numbers of stars are shown as solid circles. Two peaks are apparent, one near 
(R-I)c=0.65 and the other close to (R-I)c=1.35, the last one corresponding to the location of M35. 
Then, we have counted the number of stars in the box corresponding to the M35 peak, the box 
immediately bluer, and the two redder boxes, delineated in Figure 4 by the horizontal line segment. 
The total number of stars in all four boxes, Nj ,^(Ic), is the initial number of members of M35 in the 
considered interval of Ic. In order to remove the contamination due to field stars, we have fitted 
three Gaussian curves to Nj ,^(Ic), one corresponding to the bulk of the field stars, the other to 
M35, and the third to the transition between the former two. There is no physical reason for this 
procedure, but the shape of any particular peak is pseudo-Gaussian, and we know that Nj^^, -the 
number of stars in each box- has to be null on both sides of Figure 4 or the equivalent in other 
Ic ranges. The computed distribution of field stars is shown as a dashed line in the figure. Note 
the excellent agreement between the data and the fit for the field stars. We have subtracted the 
computed number of field stars for each box from the measured number, for the boxes corresponding 
to the location of M35. This value, N^^^^(Ic), is the number of cluster members and is listed in 
column #5 of Table 3 (in the case of the CFHT data). 

Table 3 contains the ratio between the spurious members of the cluster and the total total 
number of stars -members-|-field stars-, per bin 0.5 magnitude wide in Ic (column #4). This is 
a preliminary result, and only more accurate methods, such as proper motion measurements or 
accurate radial velocities, can guarantee the membership of a particular star or group of stars. 

4. Mass segregation in the center of the cluster 

Stellar clusters are not stable from the gravitational point of view. They disipate as they orbit 
around the Galaxy, due to the gravitational interaction with the spiral arms, interstellar clouds, 
other clusters, et cetera. In addition, since they behave as a N-body problem, they try to find 
gravitational equilibrium by expelling members, which tend to be the less massive, while the most 
massive stars move toward the center. Therefore, mass segregation appears in relaxed clusters. 

M35 has an angular diameter of, at least, 33 arcmin, since proper motion members have been 
detected at that distance (Cudworth 1971). In fact, Leonard & and Merritt (1989) stated that 
the tidal radius could reach between 33 and 66 arcmin, extending the cluster even further. Our 
survey covers an area of 27.6x27.6 arcmin, a significant fraction of the total area in the core. In 
this central part of the cluster, we have found strong evidence of a limited mass segregation, in 
a specific way. This mass segregation appears when we compare our sample with the sample of 
more massive stars from McNamara & Sekiguchi (1986a). McNamara k. Sekiguchi (1986b) found 
that essentially almost 100% of the proper motion members are within an angular radius of 20 
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arcmin, regardless of the mass (6.0<Mass<1.2 Mq), based on data from Cudworth (1971). They 
also established that there is mass segregation in this mass range, using either Cudworth's or their 
own data (sec their Figure 1). Similar conclusion has been reached by von Hippel et al. (2000), 
for stars brighter than V=19 magnitudes. However, our data, composed of less massive stars, are 
distributed homogeneously on the surveyed sky section, without any clear decrease in the number 
of stars toward the borders, indicating that the distribution of low mass stars could be extended 
far beyond the 20 arcmin radius. 

To illustrate this situation further. Figure 5 displays the cumulative distribution of stars versus 
the angular distance (in arcmin, bottom of the diagram) or radius (in parsecs, top of the diagram). 
Open circles represent the stars from Sung &; Bessell (1999) which have Ic<13 magnitudes, whereas 
solid and open triangles represent data from our survey (13<Ic<17.5 and 17.5<Ic<21.5, respec- 
tively). The mass segregation for the brighter members seems obvious. However, the two fainter 
subsamples resemble each other, except in the interval 3.5-7 arcmin, where a small difference 
appears. Moreover, the distributions in the intervals Ic=17.5-19.5 and Ic=19.5-21.5 are almost 
identical. However, some mass segregation clould be present at the very bottom of the MS, around 
the substellar limit (see Subsection 5.2.4). These facts indicate that the cluster is not fully relaxed. 



5. Luminosity and Mass Functions 

5.1. The Luminosity Function 

Figure 6 shows the Luminosity Function of M35, once the spurious members were removed. 
Solid circles represent data coming from the KPNO survey, whereas open circles correspond to the 
fainter data, collected at CFHT. We have included the errors as vertical bars (see section 5.1.1). 
Note that apparent magnitudes are indicated with the bottom of the x-axis, whereas absolute 
magnitudes appear at the top. The distribution peaks around Ic=19, as happens in the Pleiades 
(Zapatero-Osorio 1997). Fainter than this magnitude, the number of stars decreases with magni- 
tude and a strong contamination by field stars appears. Our sample is complete up to Ic=22.0 
magnitudes. Therefore, the last point is uncertain and represents only an upper limit. We have 
not considered the effect of binarity when computing the LF (or Mass Function) . 

There is an additional remark regarding the position of M35 in the CMD. We located the 
ZAMS in the diagram using the published values in the literature and appropriate conversions to 
the photometric bands wc have used. However, we did not use this information when counting the 
number of stars in M35. Therefore, our method is totally independent of distance modulus and 
reddening. However, we have obtained that the maximum density of stars, given a (R-I)c value, 
appears ^0.1 magnitudes below the ZAMS, as a visual inspection indicates in Figure 3b. This 
seems to be an effect of the low metallicity of M35. 
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5.1.1. Estimation of the errors 

To estimate the uncertainties in the Luminosity Function, we took into account the error in 
the photometry and those produced by the way we removed the contamination by field stars. The 
error in Ic is 0.05 mag. Since the bin size is 0.5 magnitudes, 20% of the number of stars -Nj_*(Ic)- 
in a particular i box could be assigned to a nearby box or i-1). In addition, wc only took into 
account the four boxes around the location of M35 (two redder and one bluer, in addition to the 
box located at the position of the cluster MS). The error in the removal of the field stars, ANgcidj 
can be estimated as the square root of the sum of the the number of stars in the external boxes 
delimiting the cluster, N™J and N^"^ in the case of the faint stars or as 30-40% for the brighter 
stars of our sample. Therefore, the total error is: 

AiV^t35(J^) = {0.22 ^ [(iV.+,^^(/^) _ AT. ,(/,))2 + (iV,,,(Ie) - iV,_i,*(Ic))'] + (AiVfield)'f ' (3) 

The values can be found in column #5 of Table 3. 

5.2. The Mass Function 

5.2.1. Estimating the stellar mass 

We estimated the masses using the Ic magnitudes. We used several isochrones from different 
groups, namely D'Antona & Mazzitelh (1994), D'Antona & Mazzitelh (1997), Baraffe et al. (1998) 
and Siess et al. (2000). In the first case, luminosities were converted into bolometric magnitudes 
and these into Ic magnitudes via bolometric corrections of Monet et al. (1992) and the distance 
modulus of the cluster. Baraffe et al. (1998) and Siess et al. (2000) provide colors and magnitudes. 
Note that the mass range of the evolutionary tracks is different for each group and that in the last 
cases they do not cover the whole range of this study. 

Figure 7a displays a comparison between the Mass Functions computed with 100 Myr isochrones 
from each model. From top to bottom, D'Antona & Mazzitclli (1994) with the magnitude-mass 
relation quoted above-, D'Antona &: Mazzitelh (1994) -where we estimated I magnitudes from 
luminosities based on Siess et al. (2000) data-, D'Antona &; Mazzitelh (1997), Baraffe et al. (1998) 
and Siess et al. (2000). For each MF, we have fitted a Salpeter law '0(M)=kxM~°, in the mass 
range 1-0.2 M©. The fits are shown as solid lines. The power law indices are, respectively, a= 0.98, 
1.10, 1.13, 0.94 and 1.14. For clarity, we have added shifts of 0.4 dex to the MFs compared to the 
previous one above it. As the figure clearly shows, the general trend is very similar in each case, 
independent of the model. Below 0.2 Mq there is a substantial decrease in the number of objects. 
Note that our completeness limit is below this value, about 0.1-0.08 M©, depending on the model. 
Therefore, this structure in the MF indicates that the efficiency of the fragmentation process of the 
original cloud of the cluster was low below 0.2 M©, producing a smaller number of objects compare 
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with stars more massive than 0.2 Mq. Another possibility is that a fraction of original low mass 
objects have been expelled during the cluster lifetime (i.e., mass segregation has occurred). 

5.2.2. The effect of the age on the Mass Function. 

As stated before, there is a considerable disagreement on the age of M35. Initial estimates were 
close to 70-100 Myr, the age of the Pleiades (Vidal 1973) and they were based on photographic 
and photoelectric photometry. More recent CCD photometry (Sung &: Bessell 1999, von Hippel et 
al. 2000) indicates that the cluster is somewhat older than the Pleiades, up to 200 Myr. Our own 
estimate, based on rotational velocity and lithium abundance of bona fide members of the cluster, 
is 175 Myr. We have tried to avoid this controversial issue by deriving the Mass Function of the 
cluster with isochrones of different ages. Figure 7b shows the results. Isochrones corresponding to 
ages of 100, 125, 160 and 200 Myr and models from Baraffc ct al. (1998) - dashed lines- and Siess 
et al. (2000) - solid lines- were used to derive the MFs. Arbitrary shifts of 0.3 and 0.5 dex were 
used to display different MFs. Again, the results are very similar in all cases. The position of the 
maximum depends slightly on the age of the isochrone and can be moved from 0.2 Mq -100 Myr- 
to 0.25 Mq -200 Myr-. In general, the effect is a shift of the Mass Function to the left side of the 
diagram, but keeping the shape of the MF unchanged. Therefore, within the uncertainty in age, 
we can ignore the effect of age on the MF. 

5.2.3. The complete Mass Function for M35 

A complete Mass Function of the open cluster M35 is displayed in Figure 8. In panels a 
and b, M35 data are shown as crosses (photometry from Sung & Bessell 1999, corrected by mass 
segregation), solid circles (Vic data from KPNO) and open circles (RIc data from CFHT). Panel a 
also includes a Salpeter law ^(M)=kxM~"^ (a=2.35, note that only stars more massive than 0.6 
Mq were used for the fit), whereas panels b and c show the M35 Mass Function in a log- log version. 
For comparison purposes, we display on this last two figures data from the Pleiades open cluster, 
as a dashed lines (Figure 8b) or thin solid lines (Figure 8c). Error bars are included for our M35 
data. 

Figure 8a clearly shows that a Salpeter law does not fit the whole Mass Function of the cluster. 
We have fitted three different power laws, yielding q;=+2.59±0.04 for the mass range 6-0.8 Mq, 
a=+0.81±0.02 for 0.8-0.2 M©, a=-0.88±0.12 for 0.2-0.08 Mq (Figure 8b). Note that the mass 
limits are different to those used in previous subsections and in Figure 7a. For the upper part of the 

main sequence of M35, Leonard & Merritt (1989) obtained a;=+2.7±0.4 (6-1 Mq). In the case of 
M35 stars slightly cooler. Sung k Bessell (1999) derived q=+2.1±0.3 (4.4-0.7 Mq). Our M35 Mass 
Function can be compared with the MF for the Pleiades, computed for mass below ~0.5 Mq, by 
Martin et al. (1998) -a=+1.00±0.15, 0.4-0.04 Mq-, Bouvier et al. (1998) -q=+0.6, 0.4-0.06 Mq-, 
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Hambly et al. (1999) -a=+0.7, 0.6-0.06 M©-, Hodgkin & Jameson (2000) -a=+0.8, 0.50-0.04 
M©-, and Martm et al. (2000) -a=+0.53, 0.4-0.06 M©-. All these MFs, arbitrarily scaled, are 
shown in Figure 8b as solid and dashed lines. In all cases, the power law indices are quite similar, 
when we compare similar mass ranges. Note that our statistics are much better due to the richness 
of the M35 cluster compared with the Pleiades. 

A comparison between a complete Mass Functions of M35 (this work) and the Pleiades (Hambly 
et al. 1999) can be found in Figure Sc. Solid circles correspond to M35 data, whereas solid triangles 
represent the Pleiades. Note the remarkable similarity between both clusters, despite the different 
age, metallicity and total mass (see next subsection) . We interpret the difference in the mass range 
0.3-0.08 Mq essentially as arising from the different binning between our data and the dataset by 
Hambly et al. (1999). Since we have a large number of candidate members at the bottom of the 
main sequence, we have been able to discriminate smaller bins close to the substellar limit, i.e., we 
have four points in the area around 0.1 Mq, where the Pleiades data presents an apparent gap, 
due to its wide bins. Therefore, we are able to distinguish more structure, and the real decrease 
of the MF in this mass range. However, we cannot rule out the possibility of a real difference in 
the MFs, even a mass segregation in the case of M35 for very low mass objects. For instance, BD 
searches in the Hyades (600-800 Myr) have not produced any positive identification (Gizis et al. 
1999). Figure 8c also includes data from Gould et al. (1997), corresponding to field stars. This 
dataset presents a flat MF in the mass range 0.5-0.08 Mq, with a step and narrow decrease at 0.1 
M0. As shown by other authors in the Pleiades and in the very young cluster (l-5Myr) associated 
to sigma Orionis (Bejar et al. 2000), the MF increases again below the substellar limit. 

In summary, it seems that the Mass Function presents four different types of behavior: From 
the top of the main sequence down to ~0.8 Mq, the MF is very step, with a power index a'^-|-2.6 
(where '0(M)=kxM~"). In the mass range 0.8-0.2 Mq, the MF is shallower, with a index close to 
a~-|-0.9. Between this point and the substellar limit, a real decreases in the M35 Mass Function 
takes place. The index can be evaluated as q;=-0.8. Finally, as shown by Bejar et al. (2000) in 
the whole substellar domain, the tendency is inverted and a new increment in the Mass Function 
appears, with a=-|-0.8. These different power law indices suggest that several mechanism (at least 
three) are acting during the fragmentation and collapse of the original cloud, competing with each 
other. 

5.2.4- The total mass of the cluster 

The total dynamical mass of M35 has been estimated as 1600-3200 Mq (Leonard & Merritt 
1989). Using mass function arguments. Sung &: Bcsscll (1999) computed a total mass of 1660 Mq. 
The average mass of a star is 0.4 Mq, which means M35 should contain from 4000 to 8000 members. 
Prior to our study, about 1500 candidate members of M35 were known. We have proposed another 
~1700, which have a contamination rate close to 30%. In total, about 2700 members of M35 have 
been identified so far. This fact seems to indicate that the total mass of the cluster is closer to the 
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1600 Mq estimate than to the 3200 Mq. In fact, our own estimates are in excellent agreement with 
this rough estimate. 

We have used the complete dataset displayed in Figure 3c to derive the total mass of the cluster. 
This sample includes our new candidate members as well as all photometric members discovered 
previously. Masses were derived using a 100 Myr isochrone (D'Antona &: Mazzitelli 1994). 

Figure 9a, 9b and 9c represent the cumulative total mass distribution against the M(Ic) absolute 
magnitude, (V-I)c,o unreddened color index and the individual mass, respectively. In the case of 
Figures 9a and c, we derived masses from the I magnitudes, whereas in the case of Figure 9b, we 
used the color. Essentially, results are very similar. The vertical long-dashed straight segment 
represents the location of the substellar limit (Figure 9b and c). The other three curves represent 
different calculations of the cumulative total mass of the cluster: The dotted line correspond to 
all probable members and candidate members. Therefore, it is affected by the field star pollution. 
The thick solid line represents the optimal estimate, where we included all probable members 
(using proper motions and radial velocities) and candidate members (our candidate members and 
previously known candidate members without known proper motions), but removing the spurious 
members (~30%). Finally, the short-dashed solid line corresponds to the probable members and 
candidate members once the field contamination has been removed. In this last case, we have not 
included candidate members from previous surveys (i.e., members with no additional information 
regarding the proper motion and/or the radial velocity) and the final value represents the minimum 
total mass of the cluster. The location of this value is indicate with a horizontal segment in Figure 
9c. 

Figure 9c clearly shows that the main contribution to the total mass of the cluster is produced 
by stars in the mass range 3-0.5 Mq, with a inflection point at 0.8 Mq. The 0.8 Mq point is also 
the location of another inflection point in the cluster MF, where the spectral index of the power 
law changes from q;=-|-2.59 to a=-|-0.81. 

The total mass within the area covered by this survey is ~2100, 1600 and 1100 Mq, for each 
case described above, where the value in the middle (1600 Mq) can be considered as a good estimate 

of the total mass of the open cluster M35 in the considered region. As a comparison, the total mass 
of the Pleiades has been estimated as 1000-1200 Mq (Meusinger et al. 1996) or 735 Mq, (Pinfield 
et al. 1998), and the total mass of Praesepe is about 626 Mq (Holland et al. 2000). Note that some 
remaining M35 members should be located in the corona of the cluster, where deep photometric 
and proper motion surveys have not been carried out. Therefore, the total mass of M35 should be 
considerably larger than the mass of the Pleiades. Pinfield et al. (1998) have estimated that the 
average mass of a Pleiades star is <mass>=0.616 Mq. In the case of M35, the average mass is 
<mass>=0.602 Mq, very similar despite the difference in the number of stars in each cluster (725 
against 2713, respectively). 



6. Summary and Conclusions 



We have obtained deep V(RI)c photometry of the open cluster M35, reaching much fainter 
stars than the previous surveys. The comparison of the location of the stars in color-magnitude 
and color-color diagrams allowed us to establish the preliminary membership of the detected stars. 
We have also estimated the contamination by field stars. Wc have analyzed this accurate, deep 
photometry of M35. Our survey has shown that in the studied color range there is a weak mass 
segregation, in opposition to the stronger segregation present for more massive stars. These data 
have allowed us to establish the Luminosity Funtion and Mass Function of the low and very low 
mass components of the clusters. The LF presents a peak at Ic=19 magnitudes, corresponding to a 
mass of ~0.3 M0. This result is very similar to the Pleiades' LF. By merging our own dataset with 
data from previous studies of the cluster, which correspond to more massive stars, wc have been 
able to obtain a complete mass function for the stellar component of the cluster, covering the mass 
range 6.0 < Mass < 0.08 Mq. When expressed in a log-log form, the MF increases monotonically in 
the 6-0.8 Mq domain with a power index of a=+2.b9. After this point, the power index is reduced 
significantly, to a value of q;=-|-0.81. Below this point, the behavior changes and the power law 
index become negative (q:=-0.8). Finally, we have estimated the total mass of the cluster in the 
area were our data were collected. After taking into account the pollution by field stars, the total 
mass contained in the central part of the cluster is '^1600 Mq. It would be very interesting to 
collect deeper photometry of the cluster, reaching well below the substellar limit, to verify whether 
the M35 MF behaves as the Pleiades and sigma Orionis Cluster or if a significant fraction of the 
M35 BDs have already been dissipated from the cluster. 
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Fig. 1. — a KPNO bright data: V magnitude against (V-I)c color, shallow exposure. Solid circles 
indicate the position of bona fide M35 members, based on radial velocity data (Barrado y Navascues 
et al. 2000a, b). The dashed and solid lines represent the completeness and detection limits, 
respectively, b KPNO bright data: V magnitude against (V-I)c color, deep exposure. Symbols as 
Figure la. c CFHT deep data: Ic magnitude against (R-I)c color index. The dashed line represent 
the completeness limit. 

Fig. 2. — Color-color magnitude diagram for M35. Solid circles represent the the M35 candidates 
selected as photometric members, whereas the crosses indicate the position of the rejected candi- 
dates (see text). As a comparison, we show IC2391 probable members (asterisks) from Barrado y 
Navascues et al. (2000c). The solid line represents the relation derived using Lcggett (1992) data, 
whereas the dashed line corresponds to the relation for M dwarfs published by Bessell (1991). 

Fig. 3. — a Final color-magnitude diagram for M35, KPNO data. Radial velocity members appear 
as solid circles (from Barrado y Navascues et al. 2000a, b), whereas other candidate members are 
displayed as open circles. Rejected candidate members, based on Figure 2, appear as crosses. 
Members from Sung & Bessell (1999) are shown as asterisk symbols, b Final M35 color-magnitude 
diagram, CFHT data. Solid circle represent stars with V(RI)c data in agreement with membership, 
whereas open circles represent candidate members which only have RI photometry, c A Complete 
color-magnitude diagram for the open cluster M35. The diagram includes all our candidate members 
and the bona fide members by Sung & Bessell (1999). Note that for the faint end (V-I)>2.0, most 
of the color indices were derived from the (R-I), using our empirical MS. 

Fig. 4. — Estimation of the pollution by field stars fitting Gaussian curves. We show one of the cuts 
along the Ic axis in the color-magnitude diagram. Three Gaussian curves were fitted as described 
in the text. The dotted line represents the fitted field star contribution. 

Fig. 5. — Mass segregation in the cluster. Open circles, solid triangles and open triangles represent 
data from Sung & Bessell (1999), KPNO and CFHT, respectively. 

Fig. 6. — Luminosity Function. Solid and open circles represent data from KPNO data and CFHT 
data, respectively. 

Fig. 7. — a Comparison of several M35 MFs derived with different models. The power law index 
for each Mass Function is shown. In all cases, 100 Myr isochrones were used. For clarity, we have 
added arbitrary shifts of 0.4 dex, respect the curve inmediately above, b Comparison of several 
M35 MFs derived with different isochrones and models. Models from Siess et al. (2000) appear 
as solid lines and circles, whereas isochrones from Baraffe et al. (1998) are plotted as dashed lines 
and open circles. For clarity, we have added arbitrary shifts of 0.3 and 0.5 dex. 

Fig. 8. — a Complete Mass Function of M35. Crosses, solid and open circles represent data from 
Leonard & Merritt (1989), KPNO data and CFHT data, respectively. The sohd line is a MF with 
a Salpeter index (a=2.35), for stars more massive than 0.6 Mq. b Complete Mass Function of 



- 21 - 



M35 in logarithmic form (thick sohd Une). The results from Leonard &; Merritt (1989) and Sung &: 
Bessell (1999) arc displayed as thin solid lines. Different Pleiades MFs, scaled with arbitrary shifs, 
arc shown as dashed lines. From top to bottom: Martin ct al. (1998), Hodgkin & Jameson (2000), 
Hambly et al. (1999), Bouvier et al. (1998) and Martin et al. (2000). The spectral power indices 
of the MF are indicated inside the parenthesis, c Complete Mass Function of M35 in logarithmic 
form (solid circles). The best fit is shown as a thick solid line. As comparison, the complete Mass 
Funtion of the Pleiades is shown as solid triangles (Hambly et al. 1999), with an arbitrary shift of 
0.4 dex. The fits and the power law indices within their range are also plotted in the Pleiades case. 
Data from Gould et al. (1997), corresponding to field stars, are included as crosses. 

Fig. 9. — a Cumulative total mass of the cluster versus the apparent Ic magnitude (top) or the 
absolute M(Ic) magnitude (bottom). Different lines represent several estimates of the total mass of 
the cluster (see text). Stellar masses were derived from the Ic magnitudes, b Cumulative total mass 
of the cluster versus the dereddened (V-I)c color index. Different lines represent several estimates 
of the total mass of the cluster (see text). The vertical long-dashed line represent the location of 
the substellar border at ~0.075 Mq for a 100 Myr old object. Stellar masses were derived from 
the (V— I)c color indices, c Cumulative total mass distribution of the cluster versus the individual 
mass. The minimum total mass for this last case is indicated with a horizontal segment. Symbols 
as in Figure 8b. Stellar masses were derived from the Ic magnitudes. 
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Table 1: Empirical M35 main sequence. We used (m-M)o=9.60, E(B-V)=0.255, E(V-I)c=0.321 
and E(R-I)c=0.178, with R v=3.0. 
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Tablc 3: Data for the LM and MF for the low mass and very low mass stars of M35 



Ic bm 


(V-I)c 


(R-l)c 


Contam. 




Mass 
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ISble 2: M35 candidates (sample). 
the authors upon request. 
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WEBDA: Open Cluster Database, Mermilliod (1996). 
MS86: McNamara & Sekiguchi (1986). 
SB99: Sung & Bessell (1999). 



